Nucleosynthesis constraint on Lorentz invariance violation in the neutrino sector 
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We investigate the nucleosynthesis constraint on Lorentz invariance violation in the neutrino sector 
which influences the formation of light elements by altering the energy density of the Universe and 
weak reaction rates prior to and during the big-bang nucleosynthesis epoch. We derive the weak 
reaction rates in the Lorentz-violating extension of the standard model. Using measurements of 
the primordial helium-4 and deuterium abundances, we give a tighter constraint on the deformed 
parameter than that derived from measurements of the cosmic microwave background anisotropics. 
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Neutrino oscillation experiments have shown that there 
arc small but non-zero mass squared differences between 
three neutrino mass cigcnstatcs (see Ref. [1] and refer- 
ence therein). However, neutrino oscillations cannot pro- 
vide absolute masses for neutrinos. Cosmology provides 
a promising way to constrain the total mass of neutrinos 
by the gravitational effect of massive neutrinos on the ex- 
pansion history near the epoch of matter-radiation equal- 
ity [2] and on the formation of large-scale structures in 
the Universe [3] (see also Ref. [4] for a review) . Recently, 
a 3(7 detection of non-zero neutrino masses is reported 
using new measurements of the cosmic microwave back- 
ground (CMB) anisotropics from the south pole telescope 
and Wilkinson microwave anisotropy probe (WMAP), 
in combination with low-redshift measurements of the 
Hubble constant, baryon acoustic oscillation feature and 
Sunyacv-Zcl'dovich selected galaxy clusters [5]. 

These observations establish the existence of physics 
beyond the standard model of particle physics. Another 
possible signal of new physics is violation of Lorentz sym- 
metry. The possibilities of Lorentz invariance violation 
were considered in quantum gravity [6] , loop gravity [7] , 
non-commutative field theory [8] , and doubly special rela- 
tivity theory [9] . Searches for Lorentz invariant violation 
with neutrinos have been performed with a wide range 
of systems [10]. Although present experiments confirm 
Lorentz invariance to a good precision, it can be broken 
in the early Universe when energies approach the Planck 
scale. Cosmological observations provide a possibility to 
test such a symmetry at high energies. 

Recently, measurements of the CMB power spectrum 
were used to probe Lorentz invariant violation in the neu- 
trino sector [11]. Lorentz invariant violation affects not 
only the evolution of the cosmological background but 
also the behavior of the neutrino perturbations. The for- 
mer alters the heights of the first and second peaks in the 
CMB power spectrum, while the latter modifies the shape 
of the CMB power spectrum. These two effects can be 
distinguished from a change in the total mass of neutrinos 
or in the effective number of neutrinos. The seven-year 



WMAP data in combination with lower-redshif measure- 
ments of the expansion rate were used to put constraints 
on the Lorentz-violating term. However, the resulting 
constraints suffer from a strong correlation between the 
Lorentz-violating term and the dark matter density pa- 
rameter [11]. 

In this letter, we use current big-bang nucleosynthe- 
sis (BBN) data to constrain Lorentz invariance violence 
in the neutrino sector. There are two effects of Lorentz 
invariant violation on BBN. The first is a correction to 
the weak reaction rate in the Lorentz-violating standard 
model extension, which governs the neutron-to-proton 
ratio at the onset of BBN. The second is a change in 
the total energy density of the Universe. Since the abun- 
dances of the light elements produced during BBN de- 
pend on the competition between the expansion rate of 
the Universe and the nuclear and weak reaction rates, the 
BBN predictions depend on the Lorentz-violating term. 
In particular the BBN-prcdictcd abundance of helium-4 
is very sensitive to the deformation parameter. 

We focus on Lorentz invariance violence only in the 
neutrino sector and consider the following deformed dis- 
persion relation 



(1) 



where E is the neutrino energy, m the neutrino mass, 
p = {p^PiY^"^ the magnitude of the 3-momentum, and 
^ the deformation parameter characterizing the size of 
Lorentz invariance violation. The dispersion relation im- 
plies that there are departures from Lorentz invariance 
in the neutrino sector if ^ 7^ 0. Such a deformed disper- 
sion relation was constructed in the framework of con- 
ventional quantum field theory [12] and derived in the 
Lorentz-violating extension of the standard model [13]. 

Here we have to point out that the dispersion rela- 
tion for neutrinos given in (1) is not very general. It ne- 
glects neutrino oscillations, possible species dependence, 
anisotropics associated with violations of rotation sym- 
metry, and CPT violation. As shown recently by Kost- 
elecky and Mewes, all of these are possible [10]. The 
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model considered in this paper is one of many possible 
Lorcntz- violating theories. 

The number density n^, and energy density pi, for mas- 
sive neutrinos with (1) are given by [11] 
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where = 2 is the number of spin degrees of freedom. 
The phase space distribution for neutrinos is the Fermi- 
Dirac distribution 



ME) = [l+cxp{E/n)] 



(4) 



where Ti, is the neutrino temperature. Thus the number 
and energy density can be written as Ji^ = {l+^)~^^'^n'i^ 
and pi, = {l+^)~'^^'^ p[^\ where ni,"' and pj?' are the stan- 
dard number and energy density, respectively. Increasing 
^ decreases both the number and energy density. The 
former leads to a reduced rate of the weak reaction prior 
to and during the BBN epoch since the reaction rate is 
proportional to the neutrino number density, while the 
latter results in a reduced expansion rate of the Uni- 
verse. Therefore, Lorentz invariant violation affects the 
nucleosynthesis of light elements. 
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FIG. 1: The ''He mass fraction (solid curve) and D/H abun- 
dance ratio (dashed curve) from BBN theory as a function 
of the deformation parameter ^ for fij/i^ = 0.0213 (the up- 
per panel) and 0.0224 (the lower panel). The vertical line 
corresponds to Lorentz invariance. 



We turn now our attention to the details of the compu- 
tation of the weak reaction rate with Lorentz invariance 
violation in the neutrino sector. At early times when the 
temperature of the Universe was T ~ 100 MeV, the num- 
ber and energy density were dominated by relativistic 
particles: electrons, positrons, neutrinos, antineutrinos 
and photons. All of the particles were kept in thermal 
equilibrium by the weak reactions 
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When the expansion rate of the Universe exceeds the 
reaction rate for n p processes, the baryons become 
uncoupled from the leptons. At the time the neutron- 
to-proton ratio is frozen, which largely determines the 
primordial helium mass fraction. To estimate the neu- 
tron abundance at the onset of BBN one has to compute 
the reaction rate. As an example, let us consider the pro- 
cess i^e + n — > ]3 + e^. The differential reaction rate per 
incident nuclcon is 
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where |A^p is the squared matrix element, to be summed 
over initial and final state spins, to„ and nip the neutron 
and proton mass, respectively, {E^, Pe) the electron four- 
momentum, and fe denotes the Fermi-Dirac statistical 
distribution for electron. The process (5) involves the 
gauge boson W as mediator. At tree level, one has 

Gp _ 

M = —j= UpJf_,{Cv - CA-f5)UnUej'^{l - Ib)^^, (9) 

where is the Fermi coupling constant, and Cy, Ca 
are the vector and axial coupling of the nucleon. After 
integrating (8) the reaction rate is 



(1- -0(1 + 0"^^ 
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(10) 



where a;*^°^ is the standard reaction rate per incident nu- 
cleon derived in [14]. The first factor on the right-hand 
side of (10) arises from the neutrino propagator in the 
Lorentz-violating extension of the standard model [15] 
and the second factor from the statistical distribution 
for neutrinos (more general Lorentz-violating corrections 
involving electrons, neutrinos, neutrons and protons were 
discussed in [16]). At tree level, the differential reaction 
rates for the other five processes in (5)-(7) can be sim- 
ply derived from (8) by properly changing the statistical 
factors and the delta function determined by the energy 
conservation for each reaction. Therefore, the corrections 
to the conversion rate of neutron into proton and its in- 
verse are the same as in (10). 
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Prom (10) we see that increasing ^ reduces the reac- 
tion rate, and thus the weak reactions freeze out at ear her 
time, corresponding to a higher freeze-out temperature. 
This leads to a larger neutron-to-baryon ratio at the onset 
of BBN and thus a larger abundance of primordial ''He 
production. On the other hand, increasing ^ also reduces 
the expansion rate of the Universe due to a decrease of 
the energy density, which means the weak reactions freeze 
out at later time without corrections to the reaction rate 
induced by the deformed parameter. This therefore re- 
sults in a lower helium-4 abundance. These two effects 
play opposite roles in the BBN prediction for the helium- 
4 abundance. The abundances of the other light nuclides 
weakly depend on ^ by changing the neutron-to-proton 
ratio and the expansion rate. 

Considering these corrections to both the reaction rate 
and the expansion rate, we now estimate the freeze- 
out temperature, Tf, determined by equating the ex- 
pansion rate with weak reaction rate. In the Friedman- 
Robertson- Walker Universe, the expansion rate obeys 
= 8nGp/3 where p oc at early times. Thus, we 
have H (X (1 — 3^/A)T^. Since the standard reaction 



rate in Eq. (10) is roughly given by uj 



(0) 



cx 



[17], 



have a; cx (1 - 15^/8)r^. Setting H ^ iiv since the free- 
ncutron decay process and its inverse are neglected at the 
BBN epoch, one derives the freeze-out temperature 
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where rj^"* is the standard one. For a large ^, the weak 
reactions freeze out at a higher temperature. This im- 
plies that effects caused by changing the reaction rate 
dominate over those by changing the expansion rate due 
to the Lorentz invariance violation in the neutrino sector. 

In order to calculate the abundances of light elements 
produced during BBN, we modified the publicly available 
PArthENoPE code [18] to appropriately incorporate the 
Lorentz-violating term in the neutrino sector. Figure 1 
shows the *He mass fraction and D/H abundance as a 
hmction of £, for flth^ = 0.0213 (the upper panel) and 
0.0224 (the lower panel). Both Yp and D/H increase as 
^ increases since the effect of changing the reaction rate 
play a leading role. Moreover, the dependence of Yp on 
^ is much larger, relative to its observational uncertain- 
ties, than that of D/H. Therefore, the primordial helium- 
4 abundance can provide a sensitive probe of neutrino 
physics with Lorentz invariance violation. 

Assuming that there are three types of neutrinos with 
vanishing chemical potentials in the Universe, the BBN- 
predicted primordial abundances depend on only two pa- 
rameters: Jlfe/i^ and ^. As shown in Figure 1, the abun- 
dance of deuterium is more sensitive to the baryon energy 
density parameter but less sensitive to the deformation 
parameter while that of helium-4 is more sensitive to ^ 
but less sensitive to Vl^h? . We use the observed primor- 
dial abundances of "^He and D in combination to constrain 
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FIG. 2: Two-dimensional joint marginalized constraints (68% 
and 95% confidence level) on the deformation parameter ^ and 
physical baryon density Q,hh? from measurements of Yp and 
D/H. The dashed line corresponds to Lorentz invariance. 



these two parameters based on the likelihood function 
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{Yp - 0.2565)2 (log[D/H] + 4.55)^ 



0.0062 



0.032 



(12) 



Here we adopt the estimate of the primordial helium 
mass fraction, Yp — 0.2565 ± 0.0060, derived in [19] using 
Monte Carlo methods to solve simultaneously for many 
possible systematic effects, based on 93 spectra of 86 low- 
metallicity extragalactic HII regions. While some have 
employed a selected subset of these data for more de- 
tailed analyses, the sources and magnitudes of system- 
atic errors have rarely been addressed. The measure- 
ment uncertainty in Yp is currently dominated by sys- 
tematic errors. For the primordial deuterium abundance, 
we use the value of log[D/H] = —4.55 ± 0.03 obtained 
in [20] from measurements of the absorption lines of seven 
quasars at high redshifts in low-metallicity hydrogen-rich 
clouds with low internal velocity dispersions. Besides 
■*He and -D, "^He and '^Li are the other two nuclides pre- 
dicted in measurable quantities by BBN. Since their post- 
BBN evolutions are complicated and their measurements 
suffer from systematical uncertainties which are difficult 
to quantify (for helium-3) or are poorly understood (for 
lithium), the observed '^He and ''Li do not provide a reli- 
able probe of BBN at present, as discussed in [21]. Thus, 
we do not include them in our constraints. 

The ''He abundance is used to provide a constrain on 
the deformation parameter while the D abundance is used 
to provide a constrain on the baryon density parameter. 
Using the combination of the ^He and D data, we find 
C = 0.034 ± 0.022 and VLbK^ = 0.0213 ± 0.0009 (68% 
confidence level). This estimated value of the deforma- 
tion parameter is consistent with Lorentz invariant ^ = 
within 95% confidence level. Compared to the results de- 
rived from the 7-year WMAP data in combination with 
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lower-rcdshift measurements of the expansion rate [11], 
BBN gives smaller uncertainties in ,f by a factor of 4 
because there is nearly no correlation between the defor- 
mation parameter and the baryon density parameter as 
shown in Figure 2. The estimate of ili,h^ is agreement 
with that from the CMB data [2] with errors. 

In summary, we have shown that the BBN puts strong 
constraint on the deformed parameter in the Lorentz- 
violating extension of the standard model, ^ = 0.034 ± 
0.022. Since the BBN-predicted abundance of helium-4 
is very sensitive to the deformed parameter but less sen- 
sitive to the baryon energy density parameter, there is 
nearly no correlation between the two parameters. Our 
results indicate no significant preference for departure 
from Lorentz symmetry in the neutrino sector in the 
early Universe. Compared to particle physics experi- 
ments, current BBN observations yield weaker constraint 
on the Lorentz-violating term. As listed in Table XIV 
of [10], previous constraints range from parts in 10^ to 
parts in 10^^ from time-of-flight measurements and vari- 
ous threshold analyses. 
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